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Abstract 

We consider the effect of disorder on coherent tunneling through two barriers 
in series, in the regime of overlapping transmission resonances. We present 
analytical calculations (using random-matrix theory) and numerical simula- 
tions (on a lattice) to show that strong mode-mixing in the inter-barrier region 
induces mesoscopic fluctuations in the conductance G of universal magnitude 
e 2 /h for a symmetric junction. For an asymmetric junction, the root-mean- 
square fluctuations depend on the ratio v of the two tunnel resistances ac- 
cording to rms G = (4e 2 //i)/3~ 1 / 2 i/(l + v)- 2 , where (3 = 1(2) in the presence 
(absence) of time-reversal symmetry. 

PACS numbers: 73.40.Gk, 72.10.Fk, 72.15.Rn 
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I. INTRODUCTION 



Resonant tunneling through two planar barriers in series is a textbook problem in quan- 
tum mechanics. Because of the separation of longitudinal and transverse motion, the prob- 
lem is essentially one-dimensional and can be solved in an elementary way. Realistic double- 
barrier junctions contain in general some amount of disorder in the region between the 
barriers. At low temperatures and small applied voltages, the inelastic electron-phonon and 
electron-electron scattering processes are suppressed, but the elastic scattering by impurities 
remains. Scattering events couple the transverse and longitudinal motion of the tunneling 
electron, which substantially complicates the problem but also leads to novel physical effects. 

The effects of disorder have been studied in the pastllHS with an emphasis on isolated 
transmission resonances (energy spacing between the resonances much greater than their 
width). Those studies are relevant for tunneling through a semiconductor quantum well, 
where the resonances are widely separated because the barrier separation L is comparable 
to the Fermi wavelength Xp. In the present paper we consider the opposite regime L ^> Ap 
of strongly overlapping resonances, relevant to metal structures (where Af is very short, 
comparable to the inter-atomic separation), or to tunneling in the plane of a two-dimensional 
electron gas (where L can be quite long, because of the large phase- coherence length). Two 
types of disorder can play a role, interface roughness at the barriers and impurities between 
the barriers. Interface roughness leads to mesoscopic (sample-to-sample) fluctuations in the 
conductance even in the absence of any phase coherence, because the tunnel probability V 
of a single barrier depends strongly on its thickness. Conductance fluctuations for a single 
rough tunnel barrier have been studied by Raikh and Ruzin.i Here we consider the case of 
impurity scattering in the absence of interface roughness. Phase coherence is then essential. 

A methodological difference with earlier work on resonant tunneling is our use of random- 
matrix theory to describe the mode-mixing in the inter-barrier region. We assume that the 
disorder is weak enough that its effect on the average conductance is negligibly small. This 
requires a mean free path I ^> YL. Still, the disorder should be sufficiently strong to fully 
mix the transverse modes in the inter-barrier region. This requires both Z <C L/r and 
W <C L/Y (where W is the transverse dimension of the junction). We may then describe 
the disorder-induced mode-mixing by a random N x N unitary matrix (N being the total 
number of propagating transverse modes at the Fermi energy). This single assumption 
permits a complete analytical solution of the statistical properties of the conductance, using 
basic results for the so-called circular ensemble of random matrices.! The circular ensemble 
is fully characterized by the symmetry index (3, which equals 1 in the presence of time- 
reversal symmetry (circular orthogonal ensemble) and 2 if time-reversal symmetry is broken 
by a magnetic field (circular unitary ensemble). (A third possibility, (5 = 4, applies to zero 
magnetic field in the presence of strong spin-orbit scattering.) 

As described in Sec. II, we find that the conductance G of the double-barrier junction 
exhibits sample-to-sample fluctuations around the classical series conductance 



(We denote by Ti and Y2 the transmission probabilities per mode through barrier 1 and 
2, and assume that these are mode-independent and < 1.) We find that the root-mean- 
square fluctuations rms G of the conductance depend only on the ratio v = ri/r 2 of the two 
transmission probabilities, according to 
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Corrections to Eq. (|1.2| ) are smaller by a factor e 2 /hG seT i es , which is 1 if iVTj ^> 1. For a 
symmetric junction [v = 1) the fluctuations are of order e 2 /h, independent of iV or r, (as 
long as NTi 3> 1). This universality is reminiscent of the universal conductance fluctuations 
in diffusive metalsfi! Just as in those systems, we expect the sample-to-sample fluctuations 
to be observable in a single sample, as reproducible fluctuations of the conductance as a 
function of Fermi energy or magnetic field. 

Eq. (|1.2|) assumes weak disorder, / TjL (but still I <C L/Ti). We generalize our 
results in Sec. Ill to stronger disorder, when the effects of the impurities on the average 
conductance have to be taken into account. As in a previous paper p where we consid- 
ered a point-contact geometry, we do this by means of the Dorokhov-Mello-Pereyra-Kumar 
(DMPK) equation.EB We find that itnpuriiy scattering leads to the appearance of a weak- 
localization effect on the average conductance (observable as a negative magnetoresistance). 
The conductance fluctuations become independent of Ti and r 2 if L ^> ICE7 1 + T^ 1 ). A 
similar conclusion was reached previously by Iida, Weidenmuller, and Zuk,H3 who studied 
the conductance fluctuations of a chain of disordered grains as a function of the coupling 
strength to two electron reservoirs. These authors found that the universal conductance 
fluctuations are recovered for a chain length L much greater than some length Lq which is 



parametrically greater than the mean free path. A more detailed comparison with Ref. |12 
is not possible, because we consider a homogeneously disordered conductor rather than a 
chain of disordered grains. 

To test our random-matrix description of mode-mixing by weak disorder, we present in 
Sec. IV results from a numerical simulation of a disordered double-barrier junction defined 
on a two-dimensional lattice. The agreement with the theory is quite reasonable. 

Two appendices to the paper contain some technical material which we need in the main 
text: In App. A we present the analogue of the Dyson-Mehta formulaEl for the circular 
ensemble, which expresses the variance of the conductance as a Fourier series. In App. B 
we discuss the application to our problem of the method of momentsBIl for the DMPK 
equation. 



II. DOUBLE-BARRIER JUNCTION WITH STRONG MODE-MIXING 

The double-barrier junction considered is shown schematically in the inset of Fig. 1. 
Since we assume Af <C L, the scattering matrix S of the whole system can be constructed 
from the scattering matrices Si of the individual barriers. The 2N x 2N unitary matrix Si 
contains two N x N submatrices fj and r\ (reflection from left to left and from right to right) 
and two other N x N submatrices tj and ^(transmission from left to right and from right 
to left). We use the polar decomposition@'0 

bi - { k < ) - \ o v t ){ r 1 / 2 -i(i - r„)v 2 ) \ o vi ) ' {2 - l) 

where the U's and Vs are N x N unitary matrices. In zero magnetic field, JJ[ = 
and V( = Vj T , so that Si is symmetric — as it should be in the presence of time-reversal 
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symmetry. The transmission matrix t of the whole system is given by 



t = t 2 {l-r' 1 r 2 Y\. 



(2.2) 



Substitution of the polar decomposition ( j2.ip yields the matrix product ttf in the form 




(2.3a) 
(2.3b) 
(2.3c) 
(2.3d) 



The eigenvalues T n of ttf are related to the eigenvalues exp(i0 n ) of Q by 



T n = (a + b cos 4> n ) 



-i 



(2.4) 



The T n 's determine the conductance G of the double-barrier junction, according to the 
Landauer formula 



where Go = 2e 2 /h is the conductance quantum. 

We consider an isotropic ensemble of double-barrier junctions, analogous to the isotropic 
ensemble of disordered wires.0 We assume that I <C L/Yi and W <C L/Yi, so that the 
tunneling is accompanied by strong mode-mixing: An electron entering the junction in mode 
n is randomly distributed among all modes m before leaving the junction. We assume in 
this section that mode-mixing is the dominant effect of the disorder, and that the reduction 
of the average conductance by the impurity scattering can be neglected. This requires 
/ ^> YiL. (The case of stronger disorder is treated in the next section.) In the polar 
decomposition ( j2.1| ) the mode-mixing is accounted for by the unitary matrices U and V. 
The number of different unitary matrices is 2f3, where (3 = 1 in zero magnetic field and 
(3 — 2 if time-reversal symmetry is broken by a magnetic field. The isotropic ensemble is 
the ensemble where the 2f3 unitary matrices are independently and uniformly distributed 
over the unitary group. In other words, the U's and V^'s are drawn independently from the 
circular unitary ensemble (CUE) of random-matrix theory! 

To determine the statistics of the conductance ( ^.5| ) we need the probability distribution 
P({0 n }) of the eigenvalues of Q. For (3 = 2, Q = U' 2 V^V[U 2 is the product of four independent 
matrices from the CUE, and hence Q is also distributed according to the CUE. For (3 = 1, 
Q = U 2 T V 1 V^ r U 2 is of the form WW T with W a member of the CUE. The ensemble of Q 
is then the circular orthogonal ensemble (COE). The distribution of the eigenvalues in the 
CUE and COE is given byi 



iv 



G — Go 22 T, 



(2.5) 



71=1 



p ({4>n}) = C Yl |exp(i0 n ) - exp(i0 m )| 



(2.6) 



where C is a normalization constant. 
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We compute the average (A) and variance Var A = (A 2 ) — (A) 2 of linear statistics 
A = Y,n=i a (0n) 011 the eigenphases 4> n . Since in the circular ensemble the 0„'s are uniformly 
distributed in (0, 2ir), the average is exactly equal to 



N r 27T 



(2.7) 



An exact expression for the variance can also be given,i but is cumbersome to evaluate. For 
N 3> 1 we can use a variation on the Dyson- Mehta formula^ (derived in App. A), 111 



Var A 



1 oo 

" y n=i 



2tt 



d0e in ^a (</>). 



(2.8a) 
(2.8b) 



For the conductance [given by Eqs. ( |2.4| ) and (|2.5|) 1, we substitute a(0) = (a + 6 cos 



with Fourier coefficients a„ = 27r(a 2 — 6 2 ) x / 2 & 



b 2)l/2 



. The results are 



(G/G ) 
Var G/G 



^(i/rx + i/r 2 - 1)- 1 
4(i-r 1 )(i-r 2 )r 2 r 



/? (r! + r 2 -rir 2 ) 



(2.9) 
(2.10) 



Equation (|2.9| ) for the average conductance is what one would expect from classical 
addition of the resistances (iVTjGo) -1 of the individual barriers. (The —1 in Eq. (|2.9|) 
corrects for a double counting of the contact resistance and becomes irrelevant for Tj 1.) 
Each member of the ensemble contains a different set of overlapping transmission resonances, 
and the ensemble average removes any trace of resonant tunneling in (G) . In a previous 
paperji^l we have shown that the average conductance differs drastically from the series 
conductance if the double-barrier junction is connected to a superconductor, but here we 
consider only normal-metal conductors. 

Eq. (|2.10|) for the conductance fluctuations tells us that VarG becomes completely in- 
dependent of iV in the limit N —>■ oo. [More precisely, corrections to Eq. ( |2.10D are of order 
(G/G )-\ which is < 1 if AT; > 1.] Since T t < 1, we may simplify Eq. (fTTDD to 



VarG/G 



r 2 r 2 

1 1 L 2 



P (ri + r 2 



(2.11) 



which depends only on the ratio rx/r 2 and not on the individual Tj's. The variance reaches 
a T-independent maximum for two equal barriers, 



VarG/Go = ±/T\ if T x = T 2 



(2.12) 



The variance is almost twice the result j=(3 1 for an isotropic ensemble of disordered 



wires,tiHi3 andnrecisely twice the result \(3 1 for an isotropic ensemble of ballistic quan- 



tum dots J 
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III. EFFECTS OF STRONG DISORDER 



In this section we relax the assumption I ^> of Sec. II, to include the case that 
the impurity scattering is sufficiently strong to affect the average conductance. We assume 
W <C L, so that we are justified in using an isotropic distribution for the scattering matrix Sl 
of the inter-barrier region.0 The scattering matrix S of the entire system is now composed 
from the three scattering matrices Si, Sl, and S2 in series. The composition is most easily 
carried out in terms of the transfer matrices M\,M L , and M 2 associated with S±,Sl, and 

respectively. The transfer matrix M of the entire system is the matrix product M = 
M 2 MiMi, so the total distribution P(M) is a convolution of the individual distributions 
Pi(M x ), P L (M L ), and P 2 (M 2 ): P = P 2 *P L * P x , where the convolution * is defined by 

P i *P j (M)= J dM' P i (MM f - 1 )P j {M f ). (3.1) 

The isotropy assumption implies that each distribution Pj(Mj) is only a function of the 
eigenvalues of MjAf/. 

We now use the fact that the convolution of isotropic distributions of transfer matrices 
commutes. (A proof is given in Ref. |9|.) This permits us to consider an equivalent system, 
with transfer matrix M = Mi l M 2 M\, where all disorder is at one side of the double-barrier 
junction — instead of in between the barriers. The L-dependence of the distribution of 
transmission eigenvalues for this system is governed by the DMPK equation, 

aS p ^ •> - g Ix, + *>W] ' (3 - 2a) 

J = njA,-A/, (3.2b) 

i<j 

where s = L/l and A n = (1 — T n )/T n . The initial condition (s — > 0) of Eq. ( p.2|) corresponds 
to taking Ml = 1, which implies for P the isotropic ensemble given by Eq. (|2.6|) . 

To compute the L-dependence of the mean and variance of the conductance, we use the 
method of moments of Mello and Stone,00 who have derived a hierarchy of differential 
equations for the moments of T q = Y^n=i ^n- The hierarchy closes order by order in an 
expansion in powers of 1/N. Mello and Stone considered a ballistic initial condition, corre- 
sponding to (Tq) N p for s — > 0. We have the different initial condition of a double-barrier 
junction. The differential equations and initial conditions for the moments are given in App. 
B. For the mean conductance and its variance we obtain 

(G/G ) = — + l(l- 2/(3) - l^ll (s 2 -s(a-p- p 2 ) + ip 3 ) , (3.3) 



s + p 3 (s + p) 3 

2 2 
15/3 + (3{s + pf 

where a has been defined in Eq. (|2.3|) and p is defined by 

P = i/r a + i/r 2 - 1. (3.5) 



Var G/G = ^ + ^-L- (s%a 2 + \p 2 - 2ap 2 + /) 

+ s{-2a 2 p + 2ap 3 - f p 5 ) + \a 2 p 2 - ip 4 - ^-p 6 ) , (3.4) 
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Corrections to Eqs. ( |3.3|) and ([3.4Q are of order (s + p)/N. For two equal barriers 
(ri = T 2 = T) in the limit r — > at fixed Ts, Eqs. ( |3.3|) and (|3.4j ) simplify to 

SG/G = (G/G„) - JV( S + p)- 1 = 1(1 - 2//3) - * ~ ^ (| + 2r s ) , (3.6) 

VarG / G »^ + M5Tr^( rV + 5 rs + l)- (37) 



Eqs. (|3.6| ) and ( |3.7|) are plotted in Fig. 1 (for (3=1). In the limit of large disorder 
(Ts ^> 1), we recover the familiar resultsEl'El for a disordered wire: 5G/G = |(1 — 2//3), 
VarG/G = (indicated by arrows in Fig. 1). In the opposite limit Ts <C 1, we find 

5G = 0, Var G/G = \p~ x — as in Sec. II [cf. Eqs. O and (gig)]. 



IV. NUMERICAL SIMULATIONS 

To test our results we have performed numerical simulations, using the recursive Green's 



function method of Ref. |22|. The disordered inter-barrier region was modeled by a tight- 
binding Hamiltonian on a two-dimensional square lattice with lattice constant d. The Fermi 
energy was chosen at 1.5 % from the band bottom, with uq = h 2 /2md 2 . Disorder was 
introduced by randomly assigning a value between ±|C/d to the on-site potential of the 
lattice points in a rectangle with L = 142d, W = 7 Id (corresponding to N = 30). We chose 
Ud = 0.6 uq, corresponding to L/l = 0.9. The transfer matrix Ml was computed numerically, 
and then multiplied with the transfer matrices Mi and Mi of the two barriers (which we 
constructed analytically, given the mode-independent tunnel probabilities Ti and T?). We 
took T 2 = 0.15 and varied Y\ between 0.05 and 0.5. These parameter values were chosen in 
order to be close to the regime TiL / L/T^ W in which disorder is expected 

to cause strong mode-mixing, without having a large effect on the average conductance (the 
regime studied in Sec. II). 

In Fig. 2 we show the comparison between theory and simulation. The solid curve is 
VarG/Go computed from 2250 realizations of the disorder potential. The dotted curve is 
the theoretical prediction from Eq. (|3.4|) for the parameter values of the simulation (and 
for (3 = 1, since there was no magnetic field). There are no adjustable parameters. The 
agreement is quite reasonable. It is likely that the remaining discrepancy is due to the 
fact that the theoretical condition iVTj ^> 1 was not well met in the simulation (where 
iVT 2 = 4.5). The value N = 30 of the simulation is already at the limit of our computational 
capabilities and we are not able to provide a more stringent numerical test of the theory. 
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APPENDIX A: DYSON-MEHTA FORMULA FOR THE CIRCULAR ENSEMBLE 

The variance Var A of a linear statistic A = J2n=i a {4>n) on the eigenphases is given by a 
double integral, 

/■2-7T /'27T" 

V&rA = - dcp d<f/a(<l>)a{<f/)K(<l),<i/), (Al) 
J o Jo 

over the two-point correlation function 

K{<t>A') = <P(0)><P(0')> - <p(0)p(0')>- (A2) 
The brackets (• • ■) denote an average over the circular ensemble, and 

N 

p(0) = £5(0-0n) (A3) 

n=l 

is the microscopic density of eigenphases. In this appendix we compute K((f>, <//) in the 
large iV-limit, using the method of functional derivatives of Ref. |23|. This leads to Eq. ( |2.8| ) 
for Var A, which is the analogue for the circular ensemble of the Dyson-Mehta formula for 
the Gaussian ensembleJi! The analogy is straightforward, but we have not found it in the 
literature.El 

We consider a generalized circular ensemble, with probability distribution 



Pv({<t>n}) =Cexp 



E^-^)+E^) > ( A4a ) 
\i<j i=i / _ 

c- 1 = r d<h r d0 2 • • • r d<p N pv^m), (a^) 

Jo Jo Jo 

U(<j>) = -ln|2sini^|. (A4c) 



The "potential" V(<p) is arbitrary If V = 0, Eq. ( |A^) is the same as the distribution ( |2.6| ) 
of the circular ensemble. The brackets (■ ■ -)y denote an average with the ^/-dependent 
distribution (|A^). Following Ref. |23], we express the two-point correlation function as a 
functional derivative of the density with respect to the potential, 

The functional derivative can be computed in the large- N limit from the relationship^ 



2t 



d(j)'U{(j) - <p')(p(<p')) v = V{4>) + const. (A6) 

Corrections to Eq. ( |A"6| ) are smaller by a factor 1/N. The additive constant is obtained from 
the normalization / d<f) (p(<j)))y = N. 

Fourier transformation of Eq. (|A6|) yields 

-^i(Pn)v = V n , n^O. (A7) 
\n\ 
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We have defined the Fourier coefficients 



fn 



2tt 



d<j>e^f(<f>), 



(A8) 



and we have used that U n = 7r/|n| for n ^ 0. From Eqs. ( |A5| ) and ( |A6| ) we see that 
K((f), (f)') = K(<ft — (f)') depends on the difference — <$' only, and is independent of V. The 
Fourier coefficients of K(<j)) are 



K n = -\n\/n(3, 



(A9a) 



for n ^ 0. Since i^o — by definition, Eq. ( |A9aj ) holds in fact for all n. Inversion of the 
Fourier transform yields the correlation function 



K(<j>) 



1 d 2 



7T 2 f3 d(f) 2 



hi 



sin 



(A9b) 



which has an integrable singularity at <fi = 0. For 0^0, K((j>) = [Att 2 (3 sin 2 (0/2)] _1 . 
Substitution of Eq. ( |A9|) into Eq. ( |AT|) gives the required analogue of the Dyson-Mehta 
formula for the large iV-limit of the variance of a linear statistic, 



VaiA 



TV 2 (3 JO 

oo 



In 



sm 



(A10) 



APPENDIX B: MOMENT EXPANSION OF THE DMPK EQUATION 

Mello and Stoned have derived from the DMPK equation fl3.2| ) a hierarchy of differential 
equations for the moments of T q = Y,n=i Tfl- The hierarchy closes order by order in the series 
expansion 



(T p ) = N*f pfl {s) + N^fM + N p ~ 2 f Pt2 {s) + . . . , (Bla) 

(T P T 2 ) = N p+1 g p+lfi (s) + NP 9p+ltl (s) + W-^M + ..., (Bib) 

(T P T 3 ) = N p+1 h p+1>0 {s) + NPh p+ltl {s) + N p -% +1>2 {s) + ..., (Blc) 

(T P T 2 2 ) = N p+2 l p+2fi (s) + NP+H p+2il (s) + N p l p+2 , 2 (s) + . . . , (Bid) 

where we have defined T = T\. For a calculation of VarG we need to determine (T p ) down 
to 0(N p ~ 2 ), (T P T 2 ) down to 0(N P ), and (T P T 3 ) and (T P T 2 2 ) only to the highest occurring 
order. The resulting set of differential equations we have to solve isEj 

f' p , (s)+pf p+lfi (s)=0, (B2a) 

g' P As) + (P + 3)^+i,o(s) = 2/ p+1|Q ( S ), (B2b) 

+Pfp+iM = (1 - 2//3) [%, (s)+pg Pt0 (sj\ , (B2c) 

£,o(s) + (P + 6)W,o(s) = % +1|Q (s), (B2d) 
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Kfi( s ) + (p + 5) Vfi.oO) = Qg P +i,o(s) - 3/ p+ i i0 (s), (B2e) 
9' P M + ( P + % p+ i,i(s) = 2f p+1>1 (s) - (1 - 2/(3) [-g' pfi (s) + 2g pfi (s) 

-4h Pj0 (s)-{p-l)l p>0 {s)), (B2f) 
f; >2 (s)+pf p+1 , 2 (s) = (1 - 2//?) [/^(s) +pg Ptl (s)_ 

+2P~ l p(p - 1) ^_i i0 (s) - Vi,o(s)] • (B2g) 

We need to determine the initial conditions f(0), g(0), h(0), and Z(0) from the distri- 
bution function (2.6) for the eigenphases in the circular ensemble. In the large- N limit, 
the linear statistic T q on the eigenphases has a Gaussian distribution with a width of or- 
der JV°. Therefore, if we write T q = (T q ) + 5T q , we know that (T q ) = O(N), (5T q ) = 0, 
((ST q ) 2n+1 ) = OiN- 1 ) and ((ST q ) 2n ) = O(N ). This implies that, for s -> 0, 



(T p ) = (T) p + \p(p - l)(T) p - 2 ((5T) 2 ) + 0(N P ~ 4 ) 
(T P T 2 ) = (T) p (T 2 ) + 0(N p - 1 ), 
(T p T 3 ) = (T) p (T 3 ) + 0(N p - 1 ), 
(T P T 2 ) = (T) P (T 2 ) 2 + 0(N P ). 

The average ((<5T) 2 ) is just VarG/G , which is given by Eq. (|2.10| ), 

((5T) 2 } =/?- 1 6 2 p" 4 + C , (^ 1 ). 

The other averages in Eq. ( P3|) follow from 

N r 2 ^ 



2vr Jo 

The resulting initial conditions read 



l\ r Z7T 

(T q ) — — / d0(a + &cos< 

Z7T JO 



/ Pl0 (o) = P - p , f P M = 0, UM = \r x vip - i)p- (p+2) b 2 

g p .o(0) = ap~( p+2 \ <7 Pj i(0) = 0, 

M°) = p-t p+2 \\a 2 p- 2 - I), / Pi0 (0) = aV P+4) - 



(B3a) 
(B3b) 
(B3c) 
(B3d) 



(B4) 



(B5) 



(B6a) 
(B6b) 
(B6c) 



The set of differential equations (|B^) can be solved by substitution of the following 
Ansatz for the p-dependence (adapted from Ref. f25|): 



= ( s + p)- (p+2 ' +n) [pV(^) +px(s) +^(s)}, 



(B7) 



where n = if x is /, n = 3 if x is g, and n = 6 if x is h or x is /. The mean and variance 
of the conductance, to order iV -1 , then follow from 



(G/G ) = Nf lt0 (s)+f ltl (s), 
Var G/G = N 2 [f 2fi (s) - f lfi (s) 2 ] + N[f 2>1 (s) - 2/ 1 , (s)/ 1 ,i(s)] + 
/ 2j 2( S )-2/ li0 ( S )/i, 2 ( S )-/i J i( S ) 2 . 



(B8) 
(B9) 



The results are Eqs. 



and (Pip . 
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FIGURES 



FIG. 1. Weak-localization correction 5G to the average conductance (in units of Go = 2e 2 /h) 
and root-mean-square fluctuations rmsG = (VarG) 1 / 2 , computed from Eqs. (3.6) and ( |3.7| ) for 
(3 = 1. The arrows give the limit TL/l ^> 1. The inset shows the geometry of the double-barrier 
junction (the disordered region is dotted). The curves plotted in the figure are for a symmetric 
junction, F\ = T2 = T <S 1. 



FIG. 2. Solid curve: variance of the conductance from a numerical simulation of an ensemble 
of disordered double-barrier junctions (L/W = 2, N = 30, s = 0.9), as a function of the ratio 
ri/r2, with T2 = 0.15 held constant. There is no magnetic field {(3 = 1). The dashed curve is the 
prediction from Eq. ( |3.4| ). There are no adjustable parameters. 
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